Analytic Study of Rotating Black-Hole Quasinormal Modes 
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A Bohr-Sommerfeld equation is derived for tlie highly-damped quasinormal mode frequencies 
uj{n ^ 1) of rotating black holes. It may be written as 2 J^{Pr + ipo) dr = {n + l/2)h, where pr 
is the canonical momentum conjugate to the radial coordinate r along a null geodesic of energy hu) 
and angular momentum hm, po = 0(tJ°), and the contour C connects two complex turning points 
of Pr- The solutions are uj{n) — —rriu) — i{(f) + n5), where {a), 5} > are functions of the black-hole 
parameters alone. Some physical implications are discussed. 
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Quantizing black holes may become an important step 
towards quantum gravity, analogous to the role played 
by atomic models in the development of quantum me- 
chanics. Thus, the "no-hair" conjecture [l| suggests that 
in a quantum theory of gravity, a black hole may be de- 
scribed by few quantum numbers related to its mass M , 
electric charge Q, and angular momentum J. The ex- 
istence of classically reversible changes in the state of a 
nonextremal black hole ^ suggests that its area A is an 
adiabatic invariant, possibly corresponding to a quantum 
entity with a discrete spectrum 

Classical black holes, like most systems with radiative 
boundary conditions, are characterized by a discrete set 
of complex ringing frequencies uj{n) — LUji + iuJi known as 
quasinormal modes (QNMs) 3]. In the spirit of Bohr's 
correspondence principle, the classical QNM spectrum 
of a black-hole should be reproduced as resonances in a 
quantum theory of gravity. QNM spectroscopy may thus 
provide valuable clues towards such a theory. In particu- 
lar, the asymptotically damped frequency ujn = ujfi{n — > 
cx)), which for a spherically-symmetric black hole de- 
pends only on the black hole parameters [e.g. may 
have a simple counterpart in quantum gravity [6|. In- 
deed, for a Schwarzschild black hole ujh — {8ttM)~^ In 3, 
such that the change in black hole entropy associated 
with AM = fuvR, AS = A{4nM'^/h) = In 3, admits a 
(triply-) degenerate quantum-state interpretation @, 0|- 
We use geometrized units where G ~ c ^ ks ^ I- 

Although uj was analytically derived for spherically 
symmetric black holes p, 01 , little is known about the 
generic and more complicated case of rotating black 
holes. Contradicting results for uj have appeared in the 
literature, although numerical convergence has recently 
been reported [8|. An analytical solution is essential in 
order to test and physically interpret these results. 



We analytically derive uj for rotating black holes in 
a method similar to the spherical black-hole analysis of 
0, by analytically continuing the relevant solution of 
Teukolsky's radial equation Q to the complex plane, and 
matching the monodromy of the wave- function along two 
different contours. Our analytical results confirm and 
generalize the numerical results of as well as admit 
a physical interpretation. In this Rapid Communication 
we outline the derivation and present the main results, 
deferring a more elaborate description of the analysis to 
a future, detailed paper. 

Teukolsky's equation. — Linear, massless field pertur- 
bations of a neutral, rotating black hole are described by 
Teukolsky's equation. For a scalar field, this equation 
can be generalized to accommodate electrically charged 
black holes O; in what follows, Q ^ is under- 
stood to apply only to such fields. The wave-function is 
separated into two ordinary differential equations using 
tP{x) = e^(™'^-'^*)S',™(cos6')i?,™(r), where x = (i,r,6',0) 
are Boyer-Lindquist coordinates. This yields radial and 
angular equations coupled by a separation constant Ai 
where Aim{ijJi — > — oo) = iAiauj + (Aq + m^) + 0(|w^ 
with yli G K [1, El- The radial equation then becomes 
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d'^ go(r)tj^ + qi{r)uj 

Q^2 + 



where Rim = A^'+^y^Ri^, A 
a = J/M, and we have defined 



0, (1) 



- 2Mr 



qa = (r' + a'Y - a' A 



(2) 



qi = -2am{2Mr - Q^) - laA^A 

+2is[r(A -I- Q2) _ M(r2 - a^)\ , (3) 



and 
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92 = -m2(A-a2)-A(s + Ao)+M2-a2-Q2 

-s(M-r)[2mm-f s(M-r)] . (4) 
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The spin-weight parameter s specifies the equation to 
gravitational (s = —2), electromagnetic (s = —1), scalar 
(s = 0), or two-component neutrino (s — —1/2) fields. 
For physical boundary conditions of purely outgoing 
waves at both spatial infinity and the event horizon (i.e. 
crossing the horizon into the black hole), Eq. ^ admits 
solutions only for a discrete set of QNM frequencies uj{n), 
where u;/ < (time decay) diverges as n ^ cxd. 

Analysis. — By defining z = V{r')dr' , with V = 
A'^{qQ + w"igi)i/2, Eq. H]) becomes 



92 

dz' 



Vi 



i? = , 



(5) 



where R = V^/^R and Vi = V"/{2V^) - Z{V'f/{W^) - 
(72/(1^ A) 2. A nonconventional tortoise coordinate z was 
defined such that the effective potential Vi — 0{\uj\^). 
The boundary condition at the horizon becomes R{r — > 
r+) ~ exp {—iLoz) cc (r — 7-^)^*'^'^+ ^ where 



Lu Res (V) 



P{cu-mn)-- + 0{\io\-^) . (6) 



Here, ft = a/ (r^ + 0^) is the angular velocity of the event 
horizon, /3 = h/iinT) = {rl + a^)/{r+ - r.), T is the 
Bekenstein- Hawking temperature, r± = M ± {AP — — 
Q^y^^ are the outer and inner horizon radii, and the tilde 
in 53 is omitted unless necessary (henceforth). R{r ~ r+) 
is multivalued, such that a clockwise rotation around r+ 
multiplies i? by a factor $1 = exp {—2TTUja+). 

Let ri and r2 — be the two complex conjugate roots 
of qo{r) lying in the fourth and in the first quadrants, 
respectively. Denote ti and t2 as the turning points of 
V [defined by V{r = ti) = 0] which lie near (a factor 
~ \uj\~^ away from) ri and r2, respectively (see Figure 
[T]). The monodromy $2 of R along a clockwise contour 
C, which passes through ti and t2 and encloses r_j_, is 
used to determine uj by demanding <I>i — ^2, as in [5]. A 
reader uninterested in details of the derivation may skip 
directly to the result, Eq. 

Near the turning points, (z — Zi) oc (r — ti)^^^ , where 
Zi = z{ti). Therefore three anti-Stokes lines, defined by 
^(iivz) = 0, emanate from ti. Two anti-Stokes lines con- 
nect ti to ^2; one (denoted ^2) crosses the real axis be- 
tween r_ and r+, while the other crosses it at r > r_|_. 
The third anti-Stokes line (h) emanating from ti extends 
to Pi, where |Pi| 00 and arg(Pi) = —tt/2. A sim- 
ilar line (/a) runs from t2 to P2, with IP2I ^ 00 and 
arg(P2) = +tt/2. A Stokes line, defined by 5(zcjz) = 0, 
emanates between every two anti-Stokes lines of ti. Let 
C be the closed, clockwise contour running from Pi to P2 
along the anti-Stokes lines h, I2 and Z3, and closing back 
on Pi through the large semicircle loo, where |r| — > cxd 
and — 7r/2 < arg(r) < 7r/2. The turning points ti and t2 
are excluded from C by partially rotating around them 
counterclockwise. Figure [T] illustrates these features in 
the r-plane. 

Along anti-Stokes lines, the WKB approximation 

R{z, Zq) ~ c+ exp [-\-iLj{z — zq)] + c_ exp [—iijj{z — zq)] 
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Im(r) 



Re(r) 



FIG. 1: Illustration of anti-Stokes (solid) and Stokes (dashed) 
lines emanating from the turning points t\ and t2 (disks) in 
the complex r-plane, for a — 0.3, Q = in the highly damped 
limit. The inner and outer horizon radii (diamonds) and com- 
ponents of the contour C are also shown. Arrows along anti- 
Stokes lines denote the direction of increasing 



holds. Off the lines, this may also be written as Cdfd + 
Csfs, where fd is exponentially large (dominant) and fg 
is exponentially small (subdominant). For < 0, the 
boundary condition at spatial infinity can be analyti- 
cally continued to Pi [5] such that P(Pi) ~ exp (-t-iwz), 
i.e. {c+,c_;2o} — {1,0; zi} up to a multiplicative fac- 
tor. This remains invariant along li till the vicinity 
of ti, so we denote R{li) = {1,0; zi}. When an anti- 
Stokes line is crossed, the dominant and subdominant 
parts exchange roles. When a Stokes line is crossed 
while circling a regular turning point, Cdfd + Cgfs be- 
comes Cdfd + (cs ± 'icd)fd, where the positive (negative) 
sign corresponds to a counterclockwise (clockwise) ro- 
tation. This so-called Stokes phenomenon [l^l implies 
that after rotating around ti from li to I2, thus crossing 
two Stokes lines and the anti-Stokes line between them, 
Rih) = {0, i;zi} — {0, i exp(— ia;(5); Z2}, where 



Z2 - Zi 



Vdr 



(7) 



Similarly, after rotating from I2 to Z3, R{lz) — 
{— exp(— 2zw(5), 0; zi}. Finally, along l^o the coefficient of 
the dominant part of the solution c+ remains invariant 
till Pi. In addition to the above changes in c+, it accu- 
mulates a phase 6+^'^'^°'+ due to the (only) singularity at 
r+ enclosed by C. Thus, the total phase accumulated by 
R along C is $2 = — exp(— 2ici;(5 -I- 2'Kijja+). For ojn > 0, 
the boundary condition at spatial infinity is continued 
to P2 and the two contours are chosen counterclockwise, 
such that the resulting equation $1 = $2 is unchanged. 

The constraint $1 = <i>2 finally yields the highly- 
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damped QNM equation [1^ 



-27rcJ(74 



-2iujS-\-27TLU(Tj^ 



(8) 



Explicitly, to order 0(|w| ^) this may be written as 



47r/3 {uj - mn) - 27ris = 2«w / y dr - 7ri(2n + 1) , (9) 



Ct., 



or in a more compact form as 



2u! I V dr — 2Tr \n + - 



(10) 




a/M 



where n e Z. Here, Ct.i (Ct.o) is a complex-plane contour 
running from ti to t2, crossing the real axis in (out) of 
the event horizon, at some point r_ < r < r_|_ (r > r+). 

Before solving for ui, note that in the highly-damped 
limit the real and the imaginary contributions to the inte- 
grals of Eqs. ((7)l- (fT0l) are easily separated. For exarnple, 
the real part of Eq. ^ may be written in the form \Vt\ 



loVr dr 



(11) 



where the complex potential Vr is given by 



qquj'^ - 2am{2Mr ~ Q'^)lu - m'^{A - a^) 
A2 



(12) 

The last term (oc taken from (72) was added to Vr 
for future use and has no effect in the highly-damped 
limit. An equation analogous to Eq. (fTT|) is found for the 
imaginary part 47r/3a;/ — 2t:s. 

QNM frequencies. — In order to obtain a closed-form 
expression for expand 2iS — Aira^ = 6o + {m5m + 'isds + 
iAiSA)cJ^'^ + 0{\uj\^^). Here 



dj = 2i 



V,dr 



with Vb — 'Zo'^^A 
[r(A + g2)-M(r2-a2)]A- 



-a(2A'fr-g2)A 



-1 -1/2 
% 



-1 -1/2 
% 



(13) 



Vs 



and Va = —% ^^^a/2. 
The integration contour Cr,o runs from ri to r2, crossing 
the real axis outside the event horizon. Since r2 = ^^i, 
{5o: ^s, ^A, ^m} are all real. Analytic expressions for these 
dj functions are readily found in terms of elliptic inte- 
grals. 

With the above definitions wc finally obtain 

Lo = —muj — + nd) , (14) 

where Q — Sm/So, S = 27r/(5o, and (f> — [sSg + AiSa — 
tt)/5q. As shown in Figures [Hand [31 these analytic results 
agree with the numerical calculations of [S] . 

Eq. yields one branch of solutions LOmin) in the 
asymptotic limit. Interestingly, in the low-n regime (and 
in spherically-symmetric black holes) two branches of so- 
lutions are identified, for given field and black-hole pa- 
rameters uM- 



FIG. 2: The real part of the highly damped QNM frequency 
uj[a) = Lijii{a;m — —1) for Q = 0, according to Eq. (|14|l (line) 
and according to the numerical results of [Sjl (circles). 



The asymptotic QNMs are not continuous at a = 



0) 



oc a 



1/3 



0, whereas 



[IJ. For Q = 0, w(a 

ujR{a = 0) = (SttM)"^ In 3. Such discontinuous behavior 
sometimes occurs in the Schwarzschild liniit. for example 
in the inner structure of the black hole Note that 

the level spacing S does continuously asymptote to the 
Schwarzschild result Au — 2TiT/h 7] as {a, Q} — > 0. 

Discussion. — We have analytically studied the highly- 
damped QNM frequencies u){n) of a rotating black hole. 
A Bohr-Sommerfeld-like equation for uj was derived 
[Eqs. (©-Uni)], analytically solved [Eq. (HI])], and shown 
to agree and generalize previous numerical results Q 
(Figures [2] and [3]). 

It is instructive to quantize the linear field perturba- 
tions described by the QNM [l^. A quantum of com- 
plex energy huj{n) and angular momentum hm may thus 
be associated with the highly-damped QNM frequency 
ujmin). Multiplying Eq. PH)) by fi yields 



pdr — 



(15) 



where p — hioV. This equation strongly resembles the 
Bohr- Sommerf eld quantization rule § pdq = (n -|- l/2)h, 




a/M 



FIG. 3: Level spacing |Aaj(a)| = 5 for Q = according to 
Eq. (|14p (line) and the numerical fit in [^ (circles). 
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where p is the canonical momentum conjugate to some 
coordinate g, and the integration is carried out along 
a closed orbit. To elucidate the connection, recall that 
the covariant radial momentum pr for geodesic motion 
of a neutral, massless particle of energy E and angular 
momentum p^, is given by 

(prA)^ = [{r^ +0^)^ ~a^A]E'^ ~2a{2Mr-Q^)Ep^ 
-{A-a^)pl-QcA , (16) 

where Qc is Carter's (fourth) constant of motion p^ . 
Comparing this with Eq. (fT2|) indicates that Vr ~ Pr, 
provided that E = Huj, p^ — hm, and Qc — 0{E'^). 
Hence, up to an 0{uj^) term which leads to an imagi- 
nary offset in uj(n), the integrand in Eq. (jlSp truly is 
of the form pdq for the above QNM quantization. The 
implied physical content of Eq. suggests that the 
full QNM spectrum may be determined by a generalized 
Bohr- Sommerf eld equation, which reduces to Eq. (|15p as 
Lui — > — oo. The general form of p is not uniquely deter- 
mined by our highly-damped analysis. Up to 0(|a;|^^) 
corrections, we may write 

p = Pr + ihsVs + iHAiVa ■ (17) 

The preceding discussion implies that Eq. psp can be 
interpreted as a complex version of the Bohr-Sommerfeld 
quantization rule. This rule was used in (the old) quan- 
tum mechanics to determine the quantum-mechanically 
allowed trajectories, as well as the quantized values of 



the associated constants of motion. Realizing the full 
meaning of Eq. (|15p may well require a quantum theory 
of gravity. Conversely, this equation can possibly be used 
to constrain and shed light on the theory. 

The quantum manifestation of a QNM may be com- 
plicated. A simple example is motivated by the outgoing 
boundary conditions of the QNMs and the symmetry of 
their frequencies W-m — —t^m evident in Eq. (fT4|) . 
These suggest that a quantum pair of opposite angular 
momentum may fundamentally correspond to a QNM; a 
positive energy quantum escaping to infinity and a neg- 
ative energy quantum falling into the black hole, in re- 
semblance of Hawking's semiclassical radiation. Under 
such circumstances, a quantum process corresponding to 
a QNM changes the black-hole mass by AM — Hlur and 
its angular momentum by AJ = hm. For such small 
changes in the black- hole parameters, the corresponding 
change in its entropy, AS = T~^{AM — flAJ), is given 
directly by Eq. (jlip . which we may now write as 



hAS = AA/4 K 2i / prdr] . (18) 



This is another indication of the adiabatic invariance of 
the area/entropy [3|. 
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